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' Abstract 

Q\ | The dc Broglie-Bohm causal theory of quantum mechanics is applied to the hydrogen atom 

(N ■ in the fully spin-dependent and relativistic framework of the Dirac equation, and in the nonrcl- 

ativistic but spin-dependent framework of the Pauli equation. Eigenstates are chosen which are 
simultaneous eigenstates of the energy H, total angular momentum M, and z component of the 
qq | total angular momentum M z . We find the trajectories of the electron, and show that in these 

^\ • eigenstates, motion is circular about the z-axis, with constant angular velocity. We compute 

the rates of revolution for the ground (n = 1) state and the n = 2 states, and show that there 
is agreement in the relevant cases between the Dirac and Pauli results, and with earlier results 
| on the Schrodinger equation. 
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1 INTRODUCTION 

s : 

In Bohm's original causal interpretation of quantum mechanics [5] , the motion of a quantum particle 
is determined by its Schrodinger wave function ip, which acts as a kind of guidance wave JjjJ. If the 
^ ■ wave function is written as 

^(x,t) = J R(x,ty 5 ^, a) 

where R and S are real- valued, then the trajectory of the particle is determined by the guidance 
relation 

p = VS. (2) 
The momentum is related to the well-known Schrodinger current j as follows, 

P = — J (3) 

P 

where p = ipty = R 2 . Comprehensive discussions of the de Broglie-Bohm causal interpretation can 
be found in and [TU] . 

It is quite natural to examine the hydrogen atom, one of the simplest quantum systems, in terms 
of the de Broglie-Bohm theory. Indeed, as originally discussed in [2], the Schrodinger guidance 
relation (j2j) predicts that p = for all real eigenstates, including the ground and all higher s 
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hydrogenic states. However, as Holland pointed out, Eq. @ is valid only for spinless particles. 
For particles with spin, the condition of Lorentz covariance on the law of motion implies that the 
momentum of a particle with spin s must be given by 

p = VS + Vlog/) x s, (4) 

where p = ip^ip [Hl^J- Only in this way can the theory be embedded in a relativistic formulation. 
Indeed, in papers [Z1IH1IH], it was shown that in order for it to be consistent with Dirac theory, the 
Schrodinger equation must be regarded as describing an electron in a definite eigenstate of spin. In 
these papers, the current vector associated with @, 

j = — p VS + — Vp x s, (5) 

m m 

was referred to as the Pauli current, the nonrelativistic limit of the Dirac current. 

The guidance law @ no longer implies that p = for real eigenstates so it is natural to ask how 
it applies to the hydrogen atom. In 0, we showed that for an electron in a spin eigenstate with 
s z = ±|, the spin-dependent term in @ will be responsible for a motion in a plane perpendicular 
to the z-axis and along a contour of constant p value. For the case of an electron in a Schrodinger 
energy /angular momentum eigenstate, V'nZmj this implies circular motion about the z-axis. 

In this paper, we examine de Broglie-Bohm trajectories for an electron in a hydrogen atom 
as described by the Pauli and Dirac equations using appropriate (spin-dependent) currents. The 
electron is assumed to be in a Pauli or Dirac eigenstate of energy and total angular momentum. Note 
that, in general, this does not imply that the electron is in a spin eigenstate of known s z value. 
We also show that under appropriate nonrelativistic limits, the angular rotation rates for Dirac 
trajectories become those of Pauli trajectories. In the cases that the electron is in a spin eigenstate 
(e.g., Is, 2s, 2po), the Pauli rotation rates agree with the Schrodinger trajectories obtained in [5]. 
The result is a coherent application of de Broglie-Bohm theory to relativistic and nonrelativistic 
hydrogen atom eigenstates. 

For both the Dirac and Pauli cases, the Schrodinger guidance formula in Q) can be generalized 
using the relationship © where j is the appropriate (Dirac/Pauli) current and p = tp^tp. First, 
consider the Dirac equation, 

d 

ih—ip = (-e<p + f3E + ct-{cv + eA))ip. (6) 
at 

Here, ip = (tpi,^,^,^) is a four-component wave function, <p and A are the scalar and vector 
potentials, Eq and p are the rest mass energy and momentum operators, e is the electric charge, 
and a and (5 are the Dirac operators. In this study, A = and the current is given by 

j = aip^aip = (j x ,jy,j z ), (7) 
where the a are the 4x4 Dirac matrices, 

-Co)- < 8 > 

and the er are the 2x2 Pauli matrices. 

If the particle is in a potential such that ecj) <C mc 2 , then there exist stationary states for 
which the average velocity v is nonrelativistic, and E ~ Eq = mc 2 . In this case, the latter two 
components of the 4-component Dirac state are smaller in magnitude than the first two components 
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by a factor of v/c. The Dirac equation may then be reduced to the Pauli equation involving the 
two components ip\ and tp2'- 

= —(-ihV + eA) 2 ^j + — <T-BiP-eV%b. (9) 
at 2m 2m 

Here ip = (^1,^2) is a two-component Pauli spinor wave function. Once again, we assume that 
A = so that the associated Pauli current is given by [3] 

j=jA+ji? = ^(V' t V^-VVV t ) + ^Vx (vW)- (io) 

Note that if one assumes that the system is in an eigenstate of the spin operator, then (fTU)) and (j3J) 
together reduce to ®. 

In the case of the Schrodinger equation for hydrogen, it is usually assumed that spin interactions 
are negligible so that the wave function can be written as a product of spatial- and spin-dependent 
terms, i.e., 

1> = il>{r,t) X ,. (11) 
As is well known, the spatial hydrogenic energy Schrodinger eigenfunctions, 

^n,e,m(r, 0, <j>) = R n/ (r)Y lm (9, (f>), (12) 

solutions to the (spinless) time-independent Schrodinger equation, are also simultaneous eigenstates 
of the orbital angular momentum operator L 2 , with eigenvalues K 2 £(£ + 1), and the operator L z , 
with eigenvalues Km. 

In the Pauli equation, where spin-orbit interactions are excluded, the orbital angular momentum 
operator I? commutes with the hamiltonian. This is not the case for the Dirac equation. For 
both the Pauli and Dirac equations, however, each component of M, the total angular momentum 
operator, commutes with the Hamiltonian H, implying that M 2 commutes with H as well. For this 
reason, it is conventional to choose eigenstates of H, M 2 and M z , with eigenvalues E n , h 2 j(j + |) 
and Km, respectively. 

There is one further subtlety: Although the orbital angular momentum does not commute with 
the hamiltonian in the Dirac case, it can be shown that £ is 'almost' a good quantum number (see 
PQ). That is to say, eigenstates can be found for which 

LV = K 2 £(£ + 1)V> + K 2 w, (13) 

where the spinor w is negligible. (Its large components actually vanish.) Hence for both the Dirac 
and Pauli cases considered in this paper, eigenstates are presented in terms of quantum numbers 
n, £, j and m for purposes of comparison, even though £ is not strictly a good quantum number in 
the Dirac case. 

Finally, in the following discussions, the time-dependent phase factor e - lE ^ t / h that accompanies 
the eigenfunctions in the solution of the time-dependent Pauli and Dirac equations will be ignored 
since it contributes nothing to the associated currents. 
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2 PAULI EIGENSTATES 

In this section we examine the Pauli current (|1U|) for some hydrogen atom eigenstates. These 
eigenstates, two-component solutions to the Pauli equation, are given by 





(14) 



Here, the R n i{r) are the standard radial wave functions for the hydrogen atom and the Y ; m _j_i (#, <j)) 

are the spherical harmonics. 1 We use spherical polar coordinates in which r is the radius, <fi is the 
angle measured counterclockwise from the x-axis and 9 is the angle measured down from the z-axis. 

As mentioned earlier, the wave functions given in (|14|) are eigenfunctions of L 2 (the orbital 
angular momentum) with eigenvalue h 2 £(£ + l), M 2 (the total angular momentum) with eigenvalue 
fr 2 j(j + 1) an d M z with eigenvalue mh. In general, however, they are not eigenstates of s z , the 
projection of spin along the z-axis. 

The eigenfunctions can be classified as follows: 

• For each n value, £ can assume the values £ = 0, 1, n — 1. 

• For each £ value, m can assume the values m = — £ + 1/2, —£ + 3/2, ...,£ — 1/2. 

• For each of the above there are two possibilities, j = £ + 1/2 and j = £ — 1/2. 

This accounts for all eigenfunctions listed in (|14|) , 

From fjlOf) . the two contributions to the velocity are given by 

_jq_ h (V>tv^-y>v^t) _ h im(^ty^) 

a p 2m e i ip^ip m e tp^tp 

and 

Vfe = k = h v x s (i6) 

where 

s = (17) 

is the 'spin vector' and m e is the mass of the electron. In order to determine trajectories for the 
above hydrogen eigenstates we must compute the velocities v a and v?, for the wave functions given 
in (ITU) . 

1 We use the following convention for the relevant functions, for consistency with [Q: 

Y lm {e,4>) = -^p lm (co S ey m4> 

v27r 



21+1 (l-m)! 1 a ro/a ri^ +m 2 / ^ n 

P £ ,_ m (a;) = (-l) m P <m (a;), m < 0. 
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We first examine the velocity v a arising from the Schrodinger current j a . Writing 

the term Imj^^V'i/'} can be shown to be 

Im^VV} = ^Imj^^d^l 2 + H 2 )e r + + v* 2 ^)e e 

-f 2 / * dv i , *dv 2 x 

+ 7^Te R ^ v ^ + v ^ )e ^- (19) 

In the above, the only nonzero term comes from the component: 

vl^=i(rn± l ~)\v k \\ k = 1,2, (20) 

so that (fT3)) yields 

™ + i(S=K)>*. Pi) 



m e rsm9\ 2 |^i| +P2I 

It remains to compute jf, and the corresponding velocity v&, with reference to (|16() . To do this, 
we first find the three components of the spin vector s = ^<jip: 

^°^ = ^T^(0 2 (2ReKt; 2 },2ImK^}> 1 | 2 -| ? ; 2 | 2 ). (22) 

From the form of the wave functions (|14[). 

2Re(«t«2) = 2c 1 c 2 N 1 N 2 P™~ l2 (e)P™ +1 2(8) cos0 (23) 

and 

2Im(^ 2 ) = 2c 1 c 2 N 1 N 2 P™~^ (e)P™^ '(0) sin0, (24) 

where the Cj are given in (|14|) and the iVj are the normalization constants of the relevant spherical 
harmonics. For simplicity of notation, define 

a = c 1 N 1 P™~K b = c 2 N 2 P™ + K (25) 

If we write 

V>W = ^-jR 2 ni (r)w, (26) 
then the vector w can be expressed in Cartesian form as 

(w x ,w y ,w z ) = (2ab cos (f>, 2ab sin <f>,a 2 — b 2 ). (27) 
We may also express w in spherical polar form, i.e., 

w x = r s sin 9 S cos 4> s 

w y = r s sin 9 S sin <p s (28) 
w z = r s cos9 s , 
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where the orientation of the spin vector s = ifr'atp is given by the angles 9 S and <p s . (See also [Q, 
p. 62-63 for a brief discussion of the spin vector.) A comparison of (|28[) with (|27jl suggests that we 
might let r s = a 2 + b 2 , <p s = (ft and then compute 9 S in terms of 9 using the relations 

a 2 -b 2 . lab . . 

cos 9 S = , sm 9 S = . (29) 

or + tr or + o 

However, this is consistent with the definition of spherical coordinates only if 2ab > since 9 S 
is restricted to the interval [0,7r]. When this condition is not met, i.e., 2ab < 0, then the polar 
coordinates for w are given by r s = a 2 + b 2 , 4> s = (p + 7r and 

q 2 -6 2 . |2o6| 2afr 

cos ^ = 2 , ia ' smgs = 2 < u2 = — 2~TT2- ( 30 ) 
or + o z a z + o z a z + tr 

In either of the above cases, the spin vector s lies in a plane defined by the position vector r 
and the z axis, which is in agreement with pQ. After some manipulation, we find that 

s = s r e r + sgeo, (31) 

where 

nr a?-b 2 . . 2ab . 
Sr = scos9(- TT¥ ) + ssm9(- TT¥ ) 

a 2 -b 2 , „, 2ab , [6) 



s d = -ssin9(— — ) + scos9( 



a 2 + b 2 a 2 +b 2 

(Here e r , eg and are the spherical polar unit vectors at the position of the electron.) Evolution of 
the position coordinates as the particle follows the trajectory implies that the spin vector precesses 
about the z-axis. This was originally described by Holland |1U| . 
From the above result we find that 

v = h V x s _ h f ^ds e 9s r \ £ 
b 2m e tp^ip 2mrs\ e dr 89 / ^ 



In other words, as was the case for v a in (|21jl. the contribution to the velocity from v& is again only 
in the direction. Therefore, for all eigenstates of the form in (|14|). the motion of the electron is 
in the direction, i.e. rotational motion about the z axis. This is in qualitative agreement with 
the Schrodinger results. 

The total speed in the direction is given by v = v a + Vb, i.e., 

ft ( . 1 ,\v 2 \ 2 - \vi\ 2 .\ . ft f ds e ds r \ 



m e rsm9\ 2 \v\\ 2 + \v2\ 2 ' 2mrs\ dr 89 

In what follows it will be useful to understand the relationship between the velocities for positive 
and negative (corresponding) values of m. Recall that for j = £ + |, 



= wr^ t{r) U/^i w»,« J (35> 



2 



and 



Y tm = -=P lm e im *, P t ,- m (x) = (-l) m P em (x). (36) 
V 2ir 
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Also, from the derivation above, the spin vector s is proportional to 

w = (2ReKw 2 },2ImK?; 2 },h| 2 - \v 2 \ 2 ). (37) 

When m is replaced by —m, we have (denoting the new term with a superscript (-) and the old 
with (+)) 

vt ] = v^^Wl = v/^^(- 1 ) m+l ^ m+ r^ m ^ ) * = v 2 +) - ( 38 ) 

Similarly, 



4~ ] = -^^V^ (39) 

= -V / ^^(- 1 ) m " l ^ m -ie- 4{m ^^ = -vl {+] . (40) 



Therefore, 

and furthermore, 



27T ' 

Ul |2_ H 2 ) (-) = _ (bl |2_ H 2 ) ( +)) (41) 



(t;J«2)H = 4 +) (-wJ (+) ) = (-«i«2)W. (42) 
All three components of w change sign when m is replaced with — m (the other eigenvalues are left 
unchanged), so that the spin vector in (|26|) changes sign. Therefore, 

,(-) . A U x MH 2 -N 2 )'-> w . _„<+) 



v -"' = ^^i m + 2< (Ni' + N')<-) ^ = - Tt ' (43> 

and 

v> ; = ; = -vl ' . (44) 

Thus, both v a and v& change sign when m changes sign, so that the overall velocity simply 
changes direction. This simplifies the computation of the rates of revolution. A similar proof holds 
for the case j = £ — | . 

Before concluding this section, we mention that in their treatment of the Pauli equation using 
Euler angles, Bohm and Schiller [I] (p. 80) deduced that the electron in a hydrogen atom eigenstate 
would execute circular motion about the principal axis with constant angular velocity. However, 
no angular velocities were computed in the paper. In the next section we compute the angular 
velocities for the first few Pauli hydrogen eigenstates. 

2.1 ANGULAR VELOCITIES FOR n = 1 and n = 2 PAULI EIGENSTATES 

We have computed explicitly the rates of revolution d<j)/dt for the first few Pauli hydrogen eigen- 
states following the procedure described above. In each case, one computes the velocity v a , followed 
by the spin vector s = i/j^ip, finding s r and sg from (|32|) . Then v& is computed to give the total 
velocity v. Since, for all cases, v points in the direction, the angular velocity d<p/dt is given by 

(45) 



dt r sin 6 

The results of our computations are presented in Tabled 

The first three results presented in Tablencorrespond to wave functions that are also eigenstates 
of s z because of the special coupling of spin and orbital angular momentum vectors. As expected, 
these rates of revolution agree with those computed in [3] for, respectively, the Is, 2s and 2y»o 
Schrodinger eigenstates. However, the final two states in Tableware not spin eigenstates. As such, 
they have no analogue in the Schrodinger case so that no comparisions of rates can be made. 
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Quantum Number n, £, j, m 


Rotation rate d<p/dt 


l,0,i,±i 


± R 

m e ar 


9 (I i +1 

u, 2 , 3= 2 


^ 2m e ar ( 1— ;f- + ^) 


2 1 ^ +5 

5 5 2> 2 


2m e «r 


o I I ±1 

5 5 2 ' 2 


±-M3 ; ) 


^5 J- 5 2' 2 


_|_ h 8cos 2 0-sin 2 
2m e ra 4 cos 2 6>+sin 2 6* 



Table 1: Angular rates of revolution for Pauli eigenstates 



3 DIRAC EIGENSTATES 



We now consider the 4-component Dirac eigenstates for hydrogen. Following 1,, they are given as 
l 

25 



follows: For j = £ + 1 



and for j = £ 



2' 



+ m+ 1 



1>i = g(r)\ r 2£ + 1 2 ^, m _i(^) 



, , £-m + l 



m + | 



2/+:] / <+i,m+^' 



2£ + 3 ' +1 ' m -r 

£ -L rji -L 2 

^4 = -»/(r)W ; - 2 Y, +1 . m+ 1 (9, 0) 



, , £-m+± 



£ -i- m -i- ! 



+ m — |. 



(46) 



(47) 



^ = -wy 2f _/ wi(^) 

// _ m _ I 

V>4 = i/wy 2 ^_ 1 2 ^- 1 , m+ |(^ 0). 

The m are the usual spherical harmonics and f(r) and <7( r ) are the normalized radial Dirac 
eigenfunctions (see p. 69). Recall that even though £ is not a good quantum number, the 
eigenstates are written in terms of £ because it is 'almost' a good quantum number, and also 
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because these solutions to the Dirac equation are built from the corresponding Pauli eigenstates. 
(For a complete discussion, see pQ.) 



3.1 TRAJECTORIES FOR GENERIC DIRAC EIGENSTATES 

In this section, we show that Bohm trajectories for Dirac hydrogen share common features. First, 
the components of the Dirac current in (J7J) may be expressed in terms of the components of the 
wave function as follows, 

-j x = 2Re{^ 4 } + 2Re{v4^ 3 } 

-J y = 2Im{v4V>4} - 2Im{v4V3} ( 4 ») 

-J z = 2Re{vIV>3} " 2Re{^ 4 }. 

We now compute these components using the hydrogenic wave functions given in ()46j) and (|47|) . 
Starting with the z component, we find that in the j '• = £ + ^ case, 



4^ = -mr)9(r)f-±^f-^^ (49) 



and 



4^ = *f{r)g{r)^ ^ + 1 2 y 2£ + 3 2 Y ttmH (e,<t>yY e+ljm+ ,(e,cl>). (50) 

Both anc ^ 4^4 are imaginary since the phases of the spherical harmonics cancel, implying 

that j z = 0. This is also the case for j = £ — ^. Therefore, in all cases, motion of the electron 
is constrained to planes of constant z. While this is a simple result, it applies to all hydrogen 
eigenstates of the forms (|46|) and (|47|) and is therefore of general interest. 

We find the other components of the current in a similar fashion. For the j ; = £ + \ case, 

~J X = 2sin0/(r) 5 (r) (J £ ^^^, m _iP Wjm+ i + / ^#^, m+ l ^ + i, m _l) • 

(51) 

We define F(cos 6) to be the quantity in brackets so that 

~j x = 2sm<Pf(r)g(r)F(cos6). (52) 

Because of its similarity in form, j y also has the F(cos 8) term: 

^j y = -2 cos (j>f(r)g(r)F (cos 9). (53) 

From © and ®, the motion of the electron in a plane of constant z is given by the following 
system of DEs: 

j x 2c/ (r)g(r)F (cos 9) sin cf> 



j y — 2cf(r)g(r)F(cos 6) cos (j) 
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From the polar forms of x and y, we have 

xx + yy = -^(x 2 + y 2 ) =0. (55) 

In other words, the motion is circular about the z-axis. 

A similar proof applies to the j ' = £ — | case. For the sake of the computations in the next 
section, the components of the current for this case are j z = 0, 

1 



= — 2 sin <j)f(r)g(r)G (cos 9) 

(56) 

and 

~i y = 2cos#(r) 5 (r)G(cos0), (57) 

where G(cos#) is the term in brackets in ()56|) . Once again we find that motion about the z-axis is 
circular. 

In summary, we have shown that electron trajectories associated with Dirac hydrogen eigenstates 
are circular, as was the case for Pauli eigenstates. In the next section, we compute some rates of 
revolution for these trajectories and their nonrelativistic limits. 

Finally, note that if m changes from positive to negative, both F(cos9) and G(cos9) simply 
undergo an overall sign change. This means that the angular rotation simply changes direction, 
but maintains the same functional form for m = ±i, as was the case for the Pauli trajectories. 

3.2 ANGULAR VELOCITIES FOR n = 1 and n = 2 DIRAC EIGENSTATES 

We find the angular rate of revolution in general from l)54j) using the relation 

y = (r sin 9 cos (f>)(j>. (58) 

From (|54|) . it follows that 

# = 2c/(r) g (r)Fcos(g) 
dt tp'ipr sin cos ^> 

Although this equation is deceptively simple in appearance, the functions /(r), g(r) - and therefore 
tp^ip - are quite complicated in form. Since we already know the quantitative motion, explicit 
computations of the rates of revolution for the general case, beyond the result given in (|59|) . are 
not particularly enlightening. 

However, it is enlightening is to examine the nonrelativistic limits of (|59|) and compare the 
results to the values computed from the Pauli equation. If the de Broglie-Bohm picture is to give 
a coherent account of the hydrogen atom, these results must agree. In what follows, we examine 
the nonrelativistic limits of ()59|) for the n = 1 and n = 2 eigenstates. Note that we compute only 
the positive m value since the angular velocity simply changes sign for negative m. 

3.2.1 n=l 

In the ground state, we have 



y{ ' \aj y 2r(2 7 i + 1) Pl (6Q) 
f(r) = -J\^9 = -Sg, 
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where 



71 = Vl-a 2 , P i = 2r/a, e x = [l + , (5 = ^-^, (61) 

a is the fine structure constant, and a is the Bohr radius. The ground state wave function is given 

by i | 

V47T V47T V47T 



This gives 

p = = — (1 + 5 2 )g 2 and F(cos 0) = sin 0. (63) 

47T 47T 

Substitution into (j59"]l yields, after cancellation, 

Tt = ^TTp- (64) 

In the nonrelativistic limit, c — > 00, which implies that a = e 2 /he — > and ji — > 1. Furthermore, 
this implies that ei — > 1 and 5^0. In order to determine the behaviour of 5c, we expand ei as 

1 a 2 



so that 



1 - ei 1 2 



1 + ei 4 

Therefore <5 — > ^<x Substitution into (fo"4'|) yields 



-a as c — > 00. (66) 



# 1 e 2 

— > -ac = —, (67) 

at r rn 



which, when written in terms of the Bohr radius a, becomes 

d6 h 



dt m P ar 



(68) 



This is the angular rotation rate for the ground state Dirac wave function. It is in agreement with 
the Schrodinger rate given by Holland [TJ3, and it is also in agreement with the rate found for the 
Pauli equation in Tabled 

3.2.2 ii=2 

1. 2S 1/2 (n = 2, I = 0, j = i, m=\) 

As in the Is case, we have the wave function 

if>l = -?==, ^2 = 0, ^3 = \=if cos 8, ip A = -L=i/ sin 0e i<? \ (69) 

V47T V47T VAtt 



where the functions /(r) and g{r) are suitably modified for the n = 2 case. Again, their 
exact functional form is not relevant, but the relationship between / and g is important; 
here, rather than (|61|) we have 



P2 = ^, JV 2 = ^2(1 + 71), ^(l + G-^— ) 2 ) , S = J\^A, (70) 
A^a V 1 + 71/ Vl + e2 
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and the number A is given by 

_ (271 + lKjVg + 2) - (iV 2 + 1>2 
(2 7l + 1)AT 2 - (jV 2 + 1) P2 

After some cancellation, substitution of the wave function into (|59|) gives, as in Q64j) . 

d 7 . ~ V^TT^' 



(71) 



(72) 



Once again, we examine how the quantities in ()7(J|) behave in the nonrelativistic limit c — > oo. 
In this case, 

e2 «l-^(T^-) 2 (73) 
2 1 + 7i 



so that 



l-£2 " 



l + e 2 4 

Prom the properties 71 — > 1, hence ./V2 — > 2, the limit of A in (J7UJ) is 



(74) 



This implies that 
and (|72*|) becomes 
which can be rewritten as 



A - ^ (75) 

2-/92 



6c -* ac W^t) (76) 



# _ e 2 (4 - p 2 ) 
eft ~~ 2hr(2 - p 2 y 



(77) 



h ( 1 

1 " ~r 7 • (78) 



dt 2m e ar \ ^ — 1 



This is the angular rotation rate for the 2s Schrodinger state given in and is also in 
agreement with the Pauli result of Table 1. 

2. 2P 3/2 (n = 2, £ = 1, j = l m = §) 
In this case the wave function is given by 

^j 1 = J—gsin e ^ fa = 0, ^3 = -* A/ — cos sin 0e^, ^ = -ifJ— sin 2 #e 3i f (79) 

V 07T V 07T V 07T 

The functions / and 5, and the relationship between them, will be the same as in the above 
case, because only m has changed. However, we now have i/j^i/j = ^- sin 2 Og 2 in the limit as 
5^0. Therefore, the expression corresponding to (|64j) is 

U -2c5 m 



dt r 

Substitution of the nonrelativistic limiting expressions from the previous case yields 

d6 h 



dt 2m P ar 



(81) 



This is the angular rotation rate for the 2pi Schrodinger state given in |5] and is also in 
agreement with the Pauli result in Table 1. 
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2P 1/2 (n = 2, 1 = 1, j = ±, m = i) 

This case is similar to the previous one although the functional forms of / and g are different. 
We have 

ip! = -^=gcos9, 7p 2 = — ^sinfle** ipz = -if-^=, ^4 = (82) 

V47T V47T V47T 



and 

w = (|)i^ < 83 > 

and most of the definitions of (|7U|) remain the same. In this case, the term A is given by 

A _ (271 + 1)^2 - (jVg ~ 1)P2 , , 

(2 7l + l)(iV 2 -2)-(iV 2 -l)p 2 - 1 j 

In the non-relativistic limit, 

A - ^ (85) 

-P2 

so that 

(it m e r 2 V 2a/ v ; 

In this case, no comparison can be made with any Schrodinger state. 

2P 3/2 (n = 2, £ = 1, j = l m = |) 
This case is somewhat different. Here, the wave function is 

1 1 



ipi = g cos 9, ip2 = j=g sinOe 1 ^ 

V27T V87T 



^ 3 = -i/^=(3cos 2 #- 1), ip 4 = -if\ —sm9cos9e ilp 

V87T V 87T 



and 

= -^(4cos 2 fl + sin 2 % 2 + — ((3cos 2 #- l) 2 + 3 sin 2 # cos 2 #)/ 2 . 
8ir 8ir 

Furthermore, the function F(cos9) becomes 

FicosO) = — sin#(8cos 2 #-sin 2 fl). 
8ir 

Again, / = — Sg, with 

£3 



and we find from (|59|) that 

d(j) 25c(8 cos 2 9 - sin 2 i 



dt r(4cos 2 9 + sin 2 0) 

As C — > OO, 



4 ' 

dd> h , 8 cos 2 — sin 2 i 



dt 2m e ar 4 cos 2 + sin 2 9 
Once again, no comparison can be made with any Schrodinger state 



(87) 



]— ~, e 3 = + , 72 = V4-a 2 (90) 

! + e3 V 72 



(91) 



1 ,a,n , 



so that 5c — ► j. After substituting and rewriting a, we obtain 

). (93) 
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In each case presented above, the nonrelativistic limit of the Dirac angular velocity agrees 
with the corresponding Pauli result given in Table ^ We expect this, since the Pauli equation 
is the nonrelativistic limit of the Dirac equation. However, the results are not obvious, since the 
expressions in (J7| and lfTU)) for, respectively, the Dirac and Pauli currents are quite different. 

4 CONCLUDING REMARKS 

In this paper, we have determined the general features of de Broglie-Bohm trajectories for en- 
ergy/total angular momentum eigenstates of the Pauli and Dirac hamiltonians for hydrogen. In 
all cases, the electron, assumed to be in an eigenstate of M z , the z-component of the total angular 
momentum, M, is confined to a plane of constant z-value and executes circular motion about the 
z-axis with a constant angular velocity d(f)/dt. As well, we have outlined a procedure to compute 
these angular velocities for general eigenstates and have explicitly computed them for the n = 1 
and n = 2 Pauli and Dirac hydrogen eigenstates. 

In the cases where the Pauli eigenstates are also eigenstates of the s z operator, our results from 
the Pauli equation agree with earlier computations of the trajectories of corresponding Schrodinger 
eigenstates [3]. Furthermore, the nonrelativistic limits of the Dirac results agree with the Pauli 
results. We have therefore shown that the de Broglie-Bohm causal picture can be applied coherently 
to the hydrogen atom, moving from the Schrodinger to the Pauli and ultimately to the Dirac 
equation. 

Finally, one may well wish to consider trajectories for Pauli or Dirac wave functions other 
than those considered in this paper. For example, it may be interesting to examine trajectories 
for particular linear combinations of eigenstates. (In jS], we examined Bohm trajectories for the 
familiar Schrodinger 2p x and 2p y orbitals used in descriptions of chemical bonding. As well, we 
examined trajectories associated with a time-varying linear combination of Is and 2pQ hydrogenic 
wave functions that simulated an electronic transition induced by an oscillating electric field.) 
The method of computing Bohm trajectories outlined in Sections 2 and 3 can be extended in a 
straightforward manner to treat such linear combinations, although the computations may well 
become quite complicated. 
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